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Abstract 

O 

I The relativistic 20-component wave equation, describing particles 

with spin 1/2 and two mass states, is analyzed. The projection op- 
erators extracting states with definite energy and spin projections, 
and density matrix are obtained. The canonical quantization of the 
field with two mass states in the formalism of the first order is per- 
(3 ■ formed and the chronological pairing of the 20-component operators 

was found. 
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O 1 Introduction 
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Ph! It is known that the problem of quark and lepton generations, and their mass 

QhI spectrum has not been solved yet. The standard model (SM) of electroweak 

interactions does not explain the number of quark and lepton generations, 
mass spectrum and contains many free parameters. Barut suggested |T] , (see 
also [2], [2, 111) to describe e-, /x-leptons by the generalized Dirac equation of 
the second order. He interpreted this equation, based on the non-perturbative 
5^ I approach to quantum electrodynamics (QED), as an effective equation for 

partly "dressed" fermions. Some investigations of the generalized Dirac equa- 
tion of the second order were made in [5| We have formulated this equation 
in the form of the first order generalized Dirac equation (FOGDE) [Hj. This 
form is convenient for analyzing symmetry properties of the theory and dif- 
ferent calculations. Here we continue to explore this matrix form of the 
equation. 

The paper is organized as follows. In Sec. 2, we formulate FOGDE 
and corresponding Lagrangian. The projection operators extracting solu- 
tions with definite energy and spin projections for free particles are given in 
Sec. 3. In Sec. 4, we perform the canonical quantization of fields. The com- 
mutation relations for 20-component functions and the chronological pairing 
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of operators were obtained. We make a conclusion in Sec. 6. Appendix 
contains some useful formulas. We use notations as in JT and the system of 
units ^ = c = 1 is chosen. 

2 First Order Field Equation 

The Barut second order generalized Dirac equations PP, describing spin-1/2 
particles, which possess two mass states, can be represented in the form 

(7.5. - + V^(x) = 0. (1) 

where = d/dxy = {d/dxm,d/d{it)), ipi^x) is a Dirac spinor and the Dirac 
matrices 7^ obey the commutation relations 7^71, + -^ul^i = '^^^lu■ The m is a 
parameter with the dimension of the mass, and a is a massless parameter in 
Eq. (1). 

Masses of spin-1/2 particles are given by 

m, = ±m[^ Ya j' = ±^ |^ Ya j" 

For a given parameter a, one can chose the positive values (for a > —1/4) of 
masses mi, m2 in Eq. (2). 

According to Barut 's idea P], the states of fermions with two masses are 
interpreted as an electron and a muon. One can apply Eq. (1) for describing 
quarks or massive neutrinos. 

It was shown in jB] that if we introduce the 20-component function 

*(.) = {^.(.)} = ( ) . o) 

where ^ 

= df,ij{x) 

m 

and index A runs values A = 0,fi, Eq. (1) becomes 

K9, + m)^(x) =0. (4) 
The 20-dimensional matrices of the first order wave equation (4) are 

a, = (^-'O + a^O'^) (g)h + eO'O ® 7,, (5) 
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where I4 is the unit 4x4 matrix, and ® means the direct product of matrices. 
The elements of the entire matrix algebra e^'^ JH] obey equations as follows: 

[e ) = OmaOn,b, e ' e ' = oabs , (oj 

where A, B, M, N = 0,1,2,3,4. The first equation in (6) shows the matrix 
elements of the matrix e*^'^, and the second equation presents the prod- 
uct of two matrices. The relativistic wave equation of the first order (4) is 
convenient for different applications because the matrices of the equation, 
a^, expressed through the elements of the entire matrix algebra with simple 
properties (6). 

The Lagrangian is given by 
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where \E'(x) = \l'"^(x)?7, 1] is the Hermitianizing matrix, and \l'"^(a;) is the 
Hermitian-conjugate wave function. The Euler-Lagrange equation (4) follows 
from the Lagrangian (7). The Hermitianizing matrix rj obeys the relations 

IZI 

rjam = -oi^r]'^, = c^tv^, = V (m = 1, 2, 3). (8) 
and is given by [HI 

= fae"''"' - ae^'^ - ® 74. (9) 



Eq. (8) guarantees that the Lagrangian (7) is the real function. By varying 
the action S = J CcPx with respect to "^{x), one may obtain the second 
Euler-Lagrange equation: 

(du^{x)^ Uu - m^{x) = 0. 

This equation may also be obtained by Hermitian conjugating Eq. (4), and 
multiplying it on i], and taking into account Eqs. (8). 



3 Mass and Spin Projection Operators 

Now we consider states of particles with definite energy, po, and momentum, 
p. It should be noted that for definite mass state, one has 

Po = \/p2 + mf, or Po = y'p2 + m|. (10) 
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This means that there is an additional quantum number r, so that r = 1, 2 
for states with masses mi and m2 jl] . We omit this index r in four momentum 
p = (p, ipo) keeping this in mind. Solutions of Eq. (4) with definite energy 
and momentum in the form of plane waves are given by 

¥ij{x) = — i==f/(±p) exp{±zpx), (11) 

where V is the normalization volume. There are two solutions (11) corre- 
sponding to Eqs. (10). Replacing Eq. (11) into Eg. (4), one obtains 

{±ip + m) U{±p) = 0, (12) 

where p = a^p^. The solution \l'^^(x) with positive frequency uj = po 
describes particles and the solution '^^Zp{x) with negative frequency uo = ~po 
corresponds to antiparticles. We use here the normalization condition 

/ ^i'^p^(x)a4^±p(x)rf^x = 1, (13) 

where '^±^{x) = (^\E'±^^(a;)^ rj. It is implied the integration over the volume 
V in Eq. (13). The difference from QED is that all functions in Eqs. (11)- 
(13) are 20-component functions and t] ^ a^. In addition, Eqs. (11)-(13) are 
valid for two mass values (r = 1,2), Eq. (10). Normalization conditions for 
20-component functions U{±p) follow from Eqs. (11)-(13): 

U{±p)a^U{±p) = -2zp„ (14) 

2 

U(±p)U(±p) = (15) 
m 

For r = 1, = —ml, and for t = 2, p"^ = —m\. Eqs. (14), (15) are analogs 
of normalization conditions for bispinors in QED. But, using the definition 
of the 20-component function ^^(x), Eq. (3), we find, from Eq. (15), the 
normalization condition for bispinor m(±p) (see also jl] for comparison): 

u{±p)u{±p) = ±^^^, (16) 

where u{±p) = u~^{±p)'y4. The normalization condition (16) for bispinors is 
different from QED. For the case of one mass state, a = 0, p"^ = —m? (see Eq. 
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(1)), we come to QED, and Eq. (16) becomes the standard normalization 
condition [7p. 

Now we construct the projection matrix extracting solutions of Eq. (12) 
corresponding to definite energy and momentum. With the help of Eqs (5), 
(6) (see Appendix), we obtain the minimal equation for the matrix p = P/^a^: 

- (1 + 2a) p^p^ + aVp = 0. (17) 

Using Eq. (17) it is not difficult to verify that the matrices 

±.ip {m ^ ip) [p^ — {1 + 2a) — rn?] . 

^ ^ 2w? [(l + 2a)p2 + 2m2] ^ ^ 

obey equations 

{±ip + m) n± = 0, (19) 

= n±, n+n_ = o. (20) 

So, matrices II-i- are projection matrices and extract solutions of Eq. (12). 
For different mass states the energy of particles, po takes two values, Eq. (10). 
In (see also |6j), we have obtained the spin projection operators 

1 / . 1\ / 2 9 



p,,,, = T5^..±5j^.^-jj, (21) 

where the operator of the spin projections on the direction of the momentum, 
p, is given by 

^abcVaJbc- (22) 



2|p 



We use the notation |p| = ypf + P2 + vl-, and the generators of the Lorentz 
group in the 20-dimensional reducible representation are j8j 

j,. = j\-:!®h+h®j^j^\ (23) 

= ^M,^ _ (24) 
4^'^ = \ (7.7. - 7.7.) • (25) 



^To have the same normahzation condition as in QED for a = Q, — — m^, one should 
make the replacement rj —i], which does not change Eqs. (8) 
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The projection operator P±i/2 obeys the equation 

(jpP±i/2 = ±\p±i/2. (26) 

Eq. (26) guarantees that the operator P±i/2 extracts spin projections s — 
±1/2. The operators p and Up commute: [p, cTj,] = and, therefore, these 
operators have the common cigcnfunction in the momentum space. 

With the aid of Eqs. (19), (26), we obtain the projection operators for 
pure spin states in the form of matrix-dyads (density matrices) 

C/,(±p) • Us{±p) = iVP±i/2n±, (27) 

where N is the normahzation constant, and we imply the matrix elements of 
matrix-dyads to be: (c/,(±p) • C7,(±p))^^ = {Us{±p))a {^^('^p))b- There is 
no summation over spin indexes s in Eq. (27). The operator C/g(±p) • C/s(±p) 
extracts the solution of Eq. (12) as well as the equation: 

aM±p) = sUsi±p), (28) 

where s = ±(1/2). The density matrices (27) are 20x20 matrices describing 
the polarization of particles for pure spin states. 
For pure spin states we use the normalization: 

2 

Us{±p)Ur{±p) = TSsr^, UMUr{-p) = 0. (29) 

Eqs. (29) generalize the normalization condition (15) on the case of states 
with definite spin projections. Prom Eq. (27), we obtain the expression for 
matrix density summed over spin projections 

^ Us{±p) ■ Us{±p) = iv±n±. (30) 

s 

Here we took into consideration the relationship (see Appendix) 

p (Pl/2 ± P-1/2) (31) 

Taking the trace in both sides of Eq. (30), and using the equahty (see 
Appendix) trll± — 2, one finds 

J2Us{±p)Us{±p) = 2N±. (32) 

s 
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If we compare Eq. (32) with the expression (29) summed over two states 
s = ±1/2, we get 

2 

N± = T-, (33) 
m 

for two values of p"^: = —m\ or = —m^. 

4 Quantization of Fermionic Fields with Two 
Mass States 

Following the general prescription ^T], we obtain from Eq. (7) the momenta: 

dC i — 

We consider here the fields \E'(x), \i/(x) as independent "coordinates" . The 
Poisson bracket {■, of the "coordinate" \I'(x) and the momentum 7r(x) is 
given by the standard equation 

{^M(x,t),7r;v(y,t)r = 5A/iv5(x-y), (36) 

and indexes M, N run 20 components. In the quantum field theory, one 
should make the substitution for fermionic fields 

{^M{^,t),nN{y,t)f ^ -z{^M(x,t),7r^(y,t)}, (37) 

where the {"^m, '^n} = '^m'^n + t^n'^m is quantum anticommutator. Replac- 
ing Eqs. (34), (35) into Eqs. (36), and taking into consideration Eq. (37), 
we arrive at 

{^M(x,t), (^(y,t)«4)^} = 25M7v5(x-y). (38) 
From Eqs. (3), (9), (38), we obtain anticommutators of Dirac spinors: 

{^^(x,t),9oV^,(y,t)} = -2z^5^,5(x-y), (39) 

{V';.(x,t),V>:(y,t)} = 0, (40) 
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where c^o = djdt, ipi^x) = ■?/'+(x)74. Commutation relations (39), (40) are 
different from those of QED because the equations considered possess higher 
derivatives, Eq. (1). In standard QED, anticommutator (40) does not equal 
zero, and there is no the equation (39). The equation like (39), including 
anticommutator of the field and its time derivative, is typical for bosonic 
fields with the replacement of anticommutators on commutators because of 
different statistics. This is due to the fact that bosonic fields obey second 
order equations as well as Eq. (1). 

The density of the Hamiltonian (the energy density) may be found from 
the equation 

n = 7i{x)do^ix) + (do^ix)) W{x) - C. (41) 

Replacing Eqs. (34), (35) into Eq. (41), and taking into consideration that 
C = for fields satisfying Eq. (4), one obtains 

n = '-'^{x)a^do^{x) - (do^ix)) a^^ix). (42) 

The density of the Hamiltonian (42) coincides with fourth components, Ti = 
T44, of the the energy-momentum tensor [H]: 

T^,u = ^ (d^'^ix)) a^^(x) - ^'^{x)a^d^^{x) 

XX a 
= ^'^{^hf.du'ipix) - 2 i^ui^i^)) 7M^(a;) + ^ {df.'^ix)) d^ipix) (43) 

-77- (d,,dutp{x)) tpix) + (duip{x)) d^4j{x) - ^ip{x)d^d^ilj{x). 
The electric current density are given by j^Oj 

j^{x) = 2^(x)a/,^(x) 

(44) 

= -i'ilj{x)-fnij{x) H ilj{x)dn'ilj{x) (dnil^lx)) ip{x). 

m m ^ ^ 

This expression includes the usual Dirac current and Barut's convective 

terms. The charge density follows from Eq. (44): jo = ""^^4 = ^E'(a;)a4\E'(x). 

Therefore, the normalization condition (13) is the normalization on the charge. 

In the second quantized theory the general solution of Eq. (4) may be 

written as 
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(45) 



Conjugated function ^'(a;) reads 



aT,p,sUr,s{p) exp(ipx) + 6+p_^C/^,,(-p) exp(-ipx) 



,(-)/ 



(46) 



where a^p^s, Or.p.s ^-re the creation and annihilation operators of particles, and 
fe+p^, &T,p,s are the creation and annihilation operators of antiparticlcs. The 
quantum number r = 1, 2 corresponds to two mass states (10). Commutation 
relations are given by 

{(J'T,p,s, (^t',p',s'} — ^ss'^tt'^pp', {(J't,p,s, Ot',p',s'} — 0, {o^p,s, (^t',p',s'} ~ 0' 

{&T,p,s, &r',p',s'} = ^ss'^rr'^pp', {&T,p,s, &t',p',s' } = 0, {6^p^^, 6^, j,, ,,, } = 0, (47) 
{ar,p,s,&T',p',s'} = {ar,p,6. = {a^p_^,6r',p',s'} = {«^p,s) ^r'.p'.s' } = 0- 

With the help of Eqs. (45)-(47), and normalization condition (13), one ob- 
tains from Eq. (42) the energy of particles-antiparticles fields 

f Hd'x = E po {a^,p,sar,p,s - K,p,sK,p,) ■ W 

T,p,S 

Like QED, we find from Eqs. (45)-(47) commutation relations for fields 



i^^rMix), ^r7v(x')} = {"^rMix), ^r7v(a;') = 0, 
{"ifrMix), "ifrNix')} = KrMN{x,x'), 

KfMNix-i — K^j^j^(x, X ) + K^j^j^(x, X ), 

KtMN{x,x')=Y.U>^+}( 



[X 



K. 



(-), 



X 



M 



M 



^iVix'] 



^i;s\x'] 



N 



N 



(49) 
(50) 

(51) 



9 



From Eqs. (45), (46), one finds 

= iUr,s{p))M (f^r,.(p))^exp[ip(x - X% 



\ _ (52) 

K^MNi^, x') = J2 iUr,s{-p))M (^r,s(-p))^ expHp(a; - x')]. 

Taking into account Eqs. (30), (33), and the relation p"^ = — m^, we obtain 
from Eqs. (52) functions as follows: 

+ , ^ ^ f ip {m - ip) [p'^ + {1 + 2a) ml - rn^] ml\ 



MN 

_ ( af^d^ (m - a^d^) [(1 + 2a) ml - rri^ - (a^S^)^] m^^ 
~ \ 2m3 [2m2 - (1 + 2a) m^] 

^?2^^"P(^^")' 

, , ^ (ip{m + ip) [p'^ + (1 + 2a) ml - m^l \ / • n 

^.M^(^) = E 4p,Vm^2m^-il + 2a)m^^] I ^"^^"^^^^ 



I (53) 

MN 



/ MAT 

_ _f (m - a^g^,) [(1 + 2a) ml - m^ - {a^d^)'^] 
~ ~ V 2m3 [2m2 - (1 + 2a) m^] J 

^^2^^"P^-^^")- 
Using the singular functions [7] 

1 1 

Ao(x) =z(A+(x) - A_(x)), 
we obtain from Eqs. (51), (53), (54) 

. / Of^d^ {m - a^(9^) [(1 + 2a) ml - m'^ - {af.df,)'^] ml ^ 
KrMN[x) - -I y 2m3 [2m2-(l + 2a)m2] ^ 

xAo(x). 



(54) 

MAT 



MN 

(55) 
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There is no summation in index r in Eqs. (49)-(55). As in QED, due to the 
properties of the function Aq{x), anticommutator {"^ m{x),^ ]\f{x')} equals 
zero if the points x and x' are separated by the space-hke interval {x — x') > 0. 
For equal times, t = t', one has {\E'm(x, 0), \E'Ar(x')} = KrMN{^ — ^', 0), where 
the function Kt-mn{^ — x', 0) may be obtained from Eq. (55) with the help 
of equalities 

ao''^Ao(x)|i=o = 0, d:i^Ao{x)\t=o = 0, doAo{x)U=o = 5(x), (56) 

where n = 1,2, 3, .... It is easy to verify, using Eq. (17), that the equations 

{a^d^ + m) K;{x) = 0, {a^d^ + m) K^{x) = (57) 

are valid. The chronological pairing of the operators in our formalism are 
given by [7] 

KM{x)KN{y) = K';MN{x-y) 

(58) 

= 0{xo- yo) K^mA^ -y)-6{yo- xo) K-^.in{x - y), 

where 9{x) is the well known theta-function. With the aid of Eqs. (53), (54), 
one finds 

KM{x)XN{y) 

(59) 

_ / Of^d^ {m - a^d„) [(1 + 2a) ml - rri^ - ml 
~ \ 2m3 [2m2 - (1 + 2a) m^] 

and the function Ac (a; — y) is given by 

A,{x -y)=9{xo- yo) A+{x -y) + e{yo- xq) A_(x - y). (60) 

5 Conclusion 

In the formalism of the first order, we have obtained the projection operators 
extracting states with definite energy and spin projections of the generalized 
Dirac equation, describing particles with spin 1/2 and two mass states. The 
density matrix was found for pure spin states. The canonical quantization 
was performed and anticommutators of 20-component fields were obtained in 
this formalism. So, FOGDE is convenient for a consideration of the conserv- 
ing currents as well as for quantization of fields. The density matrix and the 
chronological pairing of the operators found allow us to calculate different 
quantum possesses in the formalism of the first order. 



^c{x - y). 



MN 
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6 Appendix 

For convenience, we write down some matrices entering the matrix density 
and singular functions. Prom Eq. (5), one finds 

P = p^a, = /(°) (8) p + /(^) (g) h, (61) 

where p = Pu^u, and 

/(o) ^ /(I) ^ (£-'0 + ae^''^) , ^ p^p.e'^'^ (62) 

With the aid of Eq. (6), we obtain matrices as follows: 

/ = (1 + a) (g) /4 + /(^^ p + a/(2) /4, (63) 

/ = (1 + 2a) p27(°) ® p + (1 + a) p^/^^^ ®h + al^'^^ ® p, (64) 

p^ = (l + 3a + a^) p^/(°) ®/4 + (1 + 2a) p^/^^^ ®p+a (1 + a) p'^I^^^^^h, (65) 
/ = (l + 4a + 3a2) p^/(°) (g) p + (l + 3a + a^) p^/^^^ (g) /4 

(66) 

+a(l + 2a) p^/(^) (g)p, 
/ = (l + 5a + 6a2 + a^) p^/(°) (g) /4 + (l + 4a + 3a2) p^/^^) p 

(67) 

+a (l + 3a + a^) p^/^^^ (g) /4. 
The matrices (62) obey the equations: 

j(0)2 ^ j(0) ^ j(0) j(l) _^ j(l) j(0) ^ /(I) ^ 7(2)2 ^ p2j{2) ^ 

(68) 

/(O) j(2) ^ j(2) j(0) ^ /(2) j(l) _^ j(l)j(2) ^ p2j(l)_ 
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